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We have analytically solved one-dimensional interacting two-component bosonic gases with spin- 
orbit (SO) coupling by the Bethe-ansatz method. Through a gauge transformation, the effect of 
SO coupling is incorporated into a spin-dependent twisted boundary condition. Our result shows 
that the SO coupling can influence the eigenenergy in a periodical pattern. The interplay between 
interaction and SO coupling may induce the energy level crossing for the ground state, which leads 
to a transition from the ferromagnetic to antiferromagnetic state. 

PACS numbers: 67.85.-d, 67.60.Bc, 03.75.Mn 



I. INTRODUCTION 

The experimental success in manipulating cold atoms 
in effective one-dimensional (ID) waveguides has deep- 
ened our understanding of the properties of the low- 
dimensional quantum gases [H-Q- Meanwhile, studies of 
synthetic gauge field in cold atom systems have also made 
great progresses: pioneering experiments of NIST group 
have generated the effective magnetic fields, electric fields 
and spin-orbit (SO) couphng in ultracold Bose gases 
and SO-coupled Fermi gases have also been realized re- 
cently d, Q- Systems with SO coupling have revealed 
intriguing phenomena in condensed matter physics, such 
as quantum spin Hall effect and topological insulators 
B 3 J where electrons play the elemental role in these 
physical systems. The realization of SO coupling in 
cold bosonic systems opens a completely new avenue for 
studying the interesting physics beyond the traditional 
condensed matter physics. 

Many theoretical researches have investigated and pre- 
dicted interesting phenomena for SO coupled bosonic sys- 
tems. For example, a single plane wave phase or a stripe 
phase in spin-1/2 Bose-Einstein condensates (BECs) with 
SO coupling has been predicted depending on the in- 
traspecies interaction larger or smaller than that of in- 
terspecies collective modes [13 and stability 
of BECs with SO coupling have been discussed, 
and interesting phases in the presence of harmonic traps 
and rotation have also been studied [T5l - [20j . As most of 
these researches are restricted on mean-field approxima- 
tion in the weakly interacting regime, there is a rather 
basic question in what extent these results are applicable 
for an arbitrary interaction. This makes the exact solu- 
tion for SO coupled cold atom system highly desirable. 

In this paper, we analytically solve the SO coupled 
spin-1/2 bosonic gases in a ring trap by the Bethe- 



ansatz (BA) method and find that the SO coupling can 
bring some new physics to the system particularly in 
the strongly interacting regime. In the absence of SO 
coupling, the spin-1/2 bosonic system is solved analyti- 
cally by the BA method and the ground state is found to 
be the ferromagnetic state [2 1 - 24,J . However, this result 
is largely changed by adding the SO coupling. Our re- 
sults based on solving BA equations indicate that the SO 
coupling influences the eigenenergy periodically. In the 
strongly interacting regime, we find that the interplay of 
the SO coupling and interaction may produce level cross- 
ing for the ground state. Consequently, the system may 
change from a ferromagnetic ground state to an antifer- 
romagnetic state. 

The present paper is organized as follows. Section II 
introduces the model. By using a rather general trans- 
formation, SO coupling effect can be transformed to the 
twisted boundary condition. The exact solution is given 
by the BA method and the energy level crossing is shown 
in the energy spectrum. Section HI discusses the BA 
equations in the strongly interacting limit. The antifer- 
romagnetic state is proved to be the ground state in the 
infinite interaction limit. A summary is given in Section 
IV. 



II. MODEL AND SOLUTION 

We consider a two-component bosonic gas confined in 
a ID ring trap in the presence of SO coupling with the 
Hamiltonian given hy H = Hq + -ffint with 



2iaa,^]-^{x), (1) 
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where '^^ — (^|, represents two internal states of 
bosonic atoms, a denotes the spin-orbit coupling strength 
and ffz is the Pauli matrix. Here, we consider a quasi-one- 
dimensional situation with the transverse motion tightly 
confined in its ground state. The ring trap enforces 
the periodic boundary condition ^^{x) = "^aix + L). 
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The interaction term is generally represented as -ffint = 
/ dx{gih^ + g2n^ + '2gi2fi^ni), where 31^2 and 312 denote 
the strengths of intraspecies and interspecies interaction, 
which are experimentally tunable. In this work, we shall 
focus on the case with spin- independent interaction, i.e., 
gi = g2 = gi2 = c, for which interaction term can be 
represented as Hint = c/ dxn^ with h — + fij^. We 

shall set ^ = 1 in the following text for convenience. 

It is difficult to solve the Hamiltonian with spin and 
momentum coupled together. Using a rather general 
transformation [25| 

(a;) = 4, (x) 6-'^"^'^% (2) 
the Hamiltonian is rewritten as follows, 

dx^i (x) (x) ~ Na\ (3) 

and the form of Hi^t = c / dxfp' is invariant with fi — 
<h''^ and "l^^ = (<l>|,'l|). The operator ^a{x), {x') 

also satisfy commutation relation (a:) , $cr' {x)] = 
d^^iS(x — x'). Meanwhile, the total momentum now is 
represented as K = / dx^l (x) (— «^ — aaz) (x) 
which is spin-dependent. The Schrodinger equation is 
= E\^) with wave function 

I*) = E / (^1-2^^) n (^^) io) ' (4) 

K •' i=l...N 

here k denotes cri, (T2..., cr at corresponding to the spin in- 
dex for different particles. By applying the wave function 
to the Schrodinger equation, we get 



{xi, a; at) = E(I>k (xi, a; at) 



(5) 



with 



It is worth mentioning that SO effect is not vanishing 
but transformed to the spin-dependent twist boundary 
condition. Explicitly, the periodic boundary condition 
{x + L) = (a;) under the transformation Eq.Q is 
changed as 

l>,(x + L) = $,(a;)e^'^^«^, (7) 

Correspondingly, the wave function now should satisfy 

{xi, ...Xj -I- L...,xm) = e*'^^--"-^(/i„ (a;i, ...Xj...,xn) . 

(8) 

The model of two-component bosons described by 
Eq.(l6l) is analytically solvable by Bethe ansatz (BA) un- 
der periodic boundary conditions plj . The eigenstates 
can be characterized by the total spin S of the system 



which varies from to N j2. The ground state of inter- 
acting two-component bosons corresponds to the ferro- 
magnetic state with S = iV/2. Now the problem of solv- 
ing one-dimensional interacting two-component bosonic 
gases with SO coupling reduces to solve the integrable 
model of (|6|) under the twist boundary condition of ([8|), 
for which we can still obtain exact solutions by the same 
method as originally developed in (26| . The system is 
solvable by the same Bethe-type wavefunction given by 



{Xq„ Xqj^_^^ ...9 {Xq2 ) 



xAiQ,P) e'^^ ''"^'^y, 



(9) 



where Q = (72, gw) and P = {pi,p2, ■■■,Pn) de- 
note permutations of 1, N, A {Q, P) are coefficients to 
be determined, 9 {x2 — xi) is the step function and {kj} 
are quasimomenta. The coefficients fulfill the following 
relations 



AiQ;...t,j...)^Y^fAiQ;...j,z...)., 



(10) 



with Yf 

J'' 



Pn. 



k,-k,+ic ' permutates ga, qb 

in A{...qa,qb---T P)- The scattering matrix remains the 
same as in the model with periodic boundary condi- 
tions. For eigenstate with total spin •S' = ^ {N — 2M) 
(0 < M < N /2) under the twist boundary condition of 
([8]), we obtain the following Bethe-ansatz equations 



^i{kj-a)L 



N 



N , 
!=1 ^ 



M 



ki + ic Yj kj — A/3 — ic' 



kl - ic n.j 



13=1 



n k - kj - ic' 
\c - ki + id 
7— 1 ^ J 



M 



i2aL 



— — ic 



n Ac - + zc 



(11) 



(12) 



with c' = c/2, {fcj} the quasimomenta and {A^} the spin 
rapidities. From the BA equations, one can observe that 
k depends on a periodically. So we can take a — aQ + 
with ao S [O) ^] and n an integer. In fact, we need only 
consider ao S [O, f-] . The quasimomentum for ao £ 
[x' T"] '^^^ deduced from ag G [O, ^] by two steps: 
first taking — ao G f-jO] with solution —kj and — A^, 
secondly taking — ao + ^ G [x'T"] ^'^^^ solution 
unchangeable. 

Taking logarithm of the above equations, we get 



{kj — a) L = 27r/j 



N 

E 



2tan"i 



M 

-E 

/9=1 



2 tan 



-1 "-J 



(13) 



^2tan-i = 2^Jc+2aL+^2tan-i h^^. 



(14) 
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FIG. 1: (Color online) The lowest energy spectrum versus c 
for systems with A'^ = 16, a — and different total spins. The 
ground state is a ferromagnetic state and the energy levels 
satisfy E{S) < E{S') for S > S' . 



FIG. 2: (Color online) The lowest energy spectrum versus c 
for systems with A'^ = 16, a = 0.95kc/N and different total 
spins. The level crossing is shown and the ground state is the 
antiferromagnetic state in the strong interaction limit. 



where Ij with j — 1,...,N denote the density quantum 
numbers and with C = 1, M the spin quantum num- 
bers. Here Ij and are integer (half-integer) depending 
on — M is odd (even). Taking the period of a into 
account, let Ij — Ij^ — n and Jq — Jq^ — 2n, the above 
equations are reduced to 



N 



(kj — ao) L = 2TTljg — 2 tan 
1=1 

M 

-1-^2 tan 

C=i 



-1 kj hi 



-1 fcj - -^c 



(15) 



N 



E2 



tan 



-1 \ ^ ^] 



2tiJq„ + 2Q:oi 



u 



-^2tan"^ 



Ac - Xp 



(16) 



/3 = 1 



The corresponding eigenenergy is given by 

E^e- Na^ (17) 



with e — fc| . The total momentum is given by 



K = K - 2Sa, 



(18) 



with 



N M 



C=i 



In the absence of the SO couphng, the model re- 
duces to the exactly solvable interacting two-component 
bosonic model |2l| - |23| with the ground state being a fer- 
romagnetic state with 5' = N/2, which is degenerate for 
5^ G [—S, S]. For a given S, we can calculate the lowest 
eigenenergy by solving the BA equations. To give a con- 
crete example, we calculate the system with = 16 and 
show the lowest energy in a given total spin space versus 
the interaction strength in Fig.l for different total spin 
S. Here the parameters c and E are in units of kc and 
Nk"^, respectively, with kc defined as kc ~ np — ttN/L. 
Apparently, we have E{S) < E{S') for 5* > 5" in the 
whole regime of interaction strength, and the ground 
state is a ferromagnetic state with maximum S. This 
is in contrast to the case of spin-1/2 fermionic gases, for 
which Licb-Mattis theorem j27|] gives the opposite law: 
E{S) < E{S') for S < S' and the ground state is an anti- 
ferromagnetic state with 5 = 0. As the interaction tends 
to the strong interaction limit, the energies for different 
S approach to the same value and become degenerate in 
the infinite interaction limit [28l - l30| . 

In the presence of the SO coupling, the ground state 
energy spectrum for systems with different S can be also 
calculated by solving the corresponding BA equations, 
which is shown in Fig. 2 as a function of the interac- 
tion strength for a fixed SO coupling a = Q.95kc/N = 
0.957r/L. In contrast to Fig.l, the energy spectrum is 
largely changed. For weak interaction, the energies still 
satisfy the rule of E{S) < E{S') for 5" > S' . However, 
with the increase of interaction, level crossing appears 
and the ferromagnetic state is no longer the ground state. 
For large enough interaction, the energy levels fulfill the 
relation of E{S) < E{S') for S < S' . That is to say, the 
ground state of the two-component bosonic system can 
change from the ferromagnetic state to the antiferromag- 
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netic state due to the interplay between the SO couphng 
and mteraction. 



III. STRONG COUPLING LIMIT 

To see more clearly how the SO coupling affects the 
ground state energy, we discuss the strong interaction 
limit, which permits us to get some analytical expres- 
sions for the energy spectrum. In the limit of L \c\ ^ 1, 

are proportional to the interaction strength c whereas 



ki remain finite 13111 



Eq.(IT5|) and Eq.(IT6l 

plified as 



Applying the Taylor expansion to 
the equations of kj and are sim- 




2tt 



N 




2NJ 



(20) 



A A - A 
2iVtan~i -f = 27rJ(;„+2aoi+Xl2tan"i ^ (21) 



13=1 



with V — V;?^ , — TT^. For the state of M = (or 

^C — J- l/4+(A^/c)^ ^ 

even M), Ij^ = - (iV - 1) /2... {N - 1) /2, from Eq.® 



27r 



2iV/ 



with fc; 



For elementary spin excitations of 



M = 1 (odd M)[22, n = -N/2,..., N/2-1, then 



27r 

T 



— ) 

2A^/ 



Since + 5 = -^Jq ' general, the equations for quasi- 
momcntum kj are given by 



0.336 



0.334 



0.332 




FIG. 3: (Color online) The lowest energy spectrum versus a 
for systems with A'^ = 16, c/fcc — 1000 and different total 
spins. Here E is in units of Nk'^ while a is in units of kc/N. 
The energy spectrum cyclically changes in a larger scope when 
a e iO,kc). 



0.334 



0.00 




0.06 



2ttL, 



K 



In the strong coupling limit c/ p oo, kj — 2TiIj^jL + 
K/N . Without SO coupling, the quasimomenta {kj} are 
independent on the total spin. From Eq. (jl9l) . it can be 
seen that K is dependent on the total spin by the item 
2005". In contrast, with SO coupling, the quasimomenta 
{kj} are shifted by the value of 2aoS/N for various total 
spin S. Substituting the above equation into e — kj, 
we get 



N{N^ - 1) 



3 L2(i 



2p 



2N )> 



— • (23) 
N ^ ' 



From Ea. (|23p . it can be seen that the ground energy 
depends on SO coupling parameter a from three aspects. 



FIG. 4: (Color online) An enlargement of the energy spec- 
trum shown in Fig. 3 for half period of a £ (0, tt/L) is shown 
(tt/L = kc/N). There is energy crossing in this scope. The 
ferromagnetic state is the ground state while a = 0. With 
the increase of a, the ground state changes to the antiferro- 
magnetic state. 



Firstly, the spin rapidity \q is dependent on a from its 
self-consistent Eq. ((2T|) . Secondly, from Eq. (fT9)) . K is pro- 
portional to ckQ, and therefore the second item /N is 
determined by ao. Thirdly, the last item —No? is a con- 
stant for fixed a and N . In the strong interaction limit 
c/ p ^ oo, as p/c ^ 0, the contribution of spin rapidity 
can be ignored. Only the item /N will affect the en- 
ergy for different S. Here K is proportional to S from 
Ea. (IT5)) . which leads to the energy e increasing with the 
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increase of S. Therefore, SO coupling favors antiferro- 
magnetic state as the ground state. 

In the case of even N and M , the quantum numbers 
J^^ = -(M+l)/2-C where C = 1...M. From Eq. , the 
quasimomenta kj — 2'Klj^/L + 2aQS/N. From Eq.p3|. 
while c/p — > oo, the energy is given by 



k?.N 



37V2 



N 



(24) 



N 

tan" 



M 



1 \ ^ 2^ _ 2aL + y 2 tan-i 



/3=1 



In the strong interaction limit, the energy is 



^2 7v(7v2_i) 



(26) 



Apparently, the energy is lowest for 5 = 0. For large 
but finite interaction, the interplay between SO coupling 
and interaction leads to the energy level crossing. In 
Fig. 3, we show the energy spectrums versus a in the 
strongly interacting case with c/kc = 1000, where the 
SO coupling parameter is taken in the regime of a/kc G 
(0,1). The cyclical change of spectrum with the increase 
of a coincides with our previous discussion. In Fig. 4, 
the scope of SO coupling parameter a is just half the 
period a/kc S (0,7r/i) of Fig .3. With the increase of 
a, there are energy level crossing for different S. When 
a = 0, the energy levels satisfy E{S) < E{S') for S < S' . 
When a is large enough, the energy levels are opposite 
to that law. After the crossing, the larger a being, the 
larger energy difference for various S will emerge. 

Before ending the paper, we would like to give a remark 
on the interacting spin-1/2 fermionic model with SO cou- 
pling, for which the SO coupling induced level crossing 
is absent. As for the ID interacting spin-1/2 fermionic 
gas [121, the ground state is an antiferromagnetic state 
{S = 0) in the absence of SO coupling. In the presence 
of SO coupling, one can still use the gauge transforma- 
tion to transform the problem into an integrable spin-1/2 
fermionic model with a spin-dependent twist boundary 
condition [2^ [2^ , and the system is determined by the 
following BA equations 



M 



L = 2ttL 



E2 



tan 



-1 



Xr 



(25) 



Similar to the bosonic case, the SO coupling still favors 
the antiferromagnetic state as the ground state. So the 
SO coupling does not lead to level crossing for spin-1/2 
fermions, i.e., the ground state will always be an antifer- 
romagnetic state. 

IV. SUMMARY 



In summary, we have analytically solved ID interacting 
spin-1/2 bosonic gases with SO coupling. As the effect 
of SO coupling can be absorbed into the twisted bound- 
ary condition, we get the exact solution to this system by 
BA method and find that the corresponding eigenenergies 
are periodically dependent on the SO coupling. The in- 
terplay between interaction and SO coupling has revealed 
the existence of energy level crossing and the ground state 
phase transition from the ferromagnetic state to antifer- 
romagnetic state. 
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